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Abstract. An analytical approach to the one-dimensional spinless Holstein 
model is proposed, which is valid at finite charge-carrier concentrations. Spectral 
functions of charge carriers are computed on the basis of self-energy calculations. 
A generalization of the Lang-Firsov canonical transformation method is shown to 
provide an interpolation scheme between the extreme weak- and strong-coupling 
cases. The transformation depends on a variationally determined parameter that 
characterizes the charge distribution across the polaron volume. The relation 
between the spectral functions of polarons and electrons, the latter corresponding 
to the photoemission spectrum, is derived. Particular attention is paid to the 
distinction between the coherent and incoherent parts of the spectra, and their 
evolution as a function of band filling and model parameters. Results are discussed 
and compared with recent numerical calculations for the many-polaron problem. 



PACS numbers: 71.27.+a, 63.20.Kr, 71.10.Fd, 71.38.-k, 71.10.-w 
1. Introduction 

Experiments on a variety of novel materials, ranging from quasi-one-dimensional 
(ID) MX solids [1], organics [2] and quasi-2D high-Tc cuprates [3] to 3D colossal- 
magnetoresistive manganites [4, 5], provide clear evidence for the existence of polaronic 
carriers, i.e., quasiparticles consisting of an electron and a surrounding lattice 
distortion. This has motivated considerable theoretical efforts to achieve a better 
understanding of strongly coupled electron-phonon systems in the framework of 
microscopic models. 

Unfortunately, even for highly simplified models, such as the spinless Holstein 
model [6] considered here, no exact analytical solutions exist, except for the Holstein 
polaron problem with a relativistic dispersion [7] or in infinite dimensions [8]. As 
a consequence, numerous numerical studies have been carried out, focussing either 
on the empty band limit (i.e., one or two electrons only; see [9-13] and references 
therein), or on the half- filled band case in one dimension, where the Peierls transition 
takes place [14]. In contrast, very little work has been done at finite carrier densities 
away from half filling [15, 16], which are, however, often realized in experiment [1-5]. 
Recently, this so-called many-polaron problem has been addressed numerically [17- 
19]. The results have led to a fairly good understanding of many aspects, but their 
interpretation is not always straight forward, which makes analytical calculations along 
these lines highly desirable. 
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In this paper, we propose an analytical approach to the ID spinless Holstein 
model, capable of describing finite charge-carrier concentrations at arbitrary model 
parameters, including the important adiabatic intermediate-coupling (IC) regime. 

First, using standard perturbation theory based on the self-energy calculation, the 
spectral functions of charge carriers will be determined in the weak-coupling (WC) 
and the strong-coupling (SC) limits at zero temperature (T = 0). In the SC regime, 
the relation between the spectral function of polarons, determining the equilibrium 
properties (in particular, the chemical potential), and the electronic spectral fimction, 
determining the photoemission spectrum, will be discussed. Special emphasis will be 
laid on the distinction between the coherent and the incoherent parts of the spectra, 
which may be calculated separately within the present approach. 

Furthermore, using a generalization of the Lang-Firsov canonical transformation 
method [20], an interpolation scheme between the extreme WC and SC cases will be 
proposed. In particular, the canonical transformation will depend on the distance 
R characterizing the charge distribution across the polaron volume. For a given 
set of model parameters and carrier concentration n, R will be determined from 
the minimum of the total energy given by the transformed Hamiltonian in the 
first, Hartree-like approximation. With R found in this way, the polaronic and 
electronic spectral functions will be calculated to study their dependence on the model 
parameters and the carrier density. The results will be discussed with regard to recent 
numerical calculations [18, 19], which have revealed a cross-over from a system with 
polaronic carriers to a rather metallic system with increasing band filling n in the 
intermediate electron-phonon coupling regime. 

The paper is organized as follows. In section 2.1 we introduce the model used, 
whereas in section 2.2 we derive the analytical results for WC, SC and IC regimes. 
Our results are discussed in section 3, and section 4 contains our conclusions. 

2. Theory 

2.1. Model 

In this paper, we are exclusively concerned with the Holstein model (HM) of spinless 
fermions, which describes electrons coupled locally to Einstein phonons. Although we 

shall use different canonical transformations to describe the IC and SC regime later 
on, the Hamiltonian can be written in the general form 



where the definition of jj and Cij will be different depending on the approach used. 

Here c| (c-) creates (annihilates) a spinless fermion at site i, &| and are bosonic 
operators for the dispersionless phonons of energy ojq {h = 1), and the strength of 
the electron-phonon interaction is specified by the dimensionless coupling constants 
A = Ep/2t and g = ^ -Ep/tco in the adiabatic {loqII <C 1) and anti-adiabatic (i^o/i ^ 1) 
regimes, respectively, where E-^ is the well-known polaron binding energy in the atomic 
limit [Cij = for z ^ j in equation (1)]. 

In the WC case, in which we use the original, untransformed Holstein 
Hamiltonian, we have r\ = where \x denotes the chemical potential, and non-zero 
coefiicients 




(1) 



(2) 
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In contrast, the starting point in the SC regime will be the Hamiltonian with 
Cu=0, Cuj)=te-s^'''-''i-'>'^+''^K (3) 



T] = —g^ujQ — 1-1= —Ep — and 



2.2. Green functions approach 

We treat the HM (1) using the formalism of the generalized Matsubara Green functions 
introduced by Kadanoff and Baym [21] and Bonch-Bruevich and Tyablikov [22], and 
applied to the single-polaron problem by Schnakenberg [23]. The Green fimction 
equation of motion deduced from H may be converted into an equation for the self- 
energy of the spinless fermions, and can be solved by iteration (see [23-25] for details). 
In the second iteration step, the self-energy is obtained in the form 

S(jlTi; j2T2) = ~{Cj,j.,)S{ti - T2) 

j'j" 

where G^{jiTi;j2T2) represents the first-order fermionic Matsubara Green function, 
and the symbol denotes the time ordering operator acting on the imaginary times . 
Fourier-transforming both sides of equation (4), carrying out the standard summation 
over the Matsubara boson frequencies [26], and using the analytical continuation in 
the complex frequency plane, the retarded momentum- and energy-dependent fermion 
Green function follows as 

where is the fermionic band dispersion in the first approximation, lu = uj + ie 
(e 0+), and ^{k, ui) denotes the Fourier transform of the coUisional part of the self- 
energy given by the second term on the r. h. s of equation (4). The related, normalized 
spectral function is given by 

A{k, co) = --ImG^{k, CO + iO+). (6) 

Although the insulating Pcicrls phase with long-range charge-density-wave order, 
which is the ground state of the half-filled spinless Holstein model above a critical 
coupling strength depending on ujq [14], can in principle be incorporated into the 
present theory, this has not been done here, as we are interested in intermediate 
band fillings away from n = 0.5. Furthermore, we neglect any extended pairing for 
< n < 0.5, as well as the possible formation of a polaronic super lattice in the half- 
filled band case (n = 0.5). Finally, the possibility of phase separation, which can in 
principle occur in the present model, is not considered. 



2.2.1. Weak coupling In the WC limit the self-energy is given by 

E(fc,aJ) 



_._LOaEp d^ [ l-nF(C-M) , »^f(^-m) 



7tw j_w VI - {^/wy 



(7) 



-UJO- H) LO + LOO - - IJ.) 

with the Fermi function n-p{Lo) = (c'^" + 1)~^ and the bare bandwidth W = 2t. At 
T = 0, and defining ^0 — max (/i, — W), we get 



ReS(fc,a;) '^"^^ 



ttW 



w 



Vl-i^/Wy^+too-iC- fx) 
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(8) 



and 

Im w) 



+ 



5[u; - Wo - (C - m)] 



(9) 



In the sequel, we shall distinguish between coherent and incoherent contributions 
to the single-particle spectral functions, as defined by a zero and non-zero imaginary 
part of the self-energy, respectively. The coherent part of the spectrum is given by 

A''{k,uj) = Zk5[uj-{Ek-ii)], (10) 

where the renormalized band energy is the solution of 



Ek = £,k + 



-kW 



V 



1 



w Vl - Ek^ujQ-i 

^ d^ 1 



4 Vr^ 



(^/T4^)2 - Wo - e 



(11) 



with = — M^cosfc (the bare band dispersion), and the spectral weight takes the 
form 



1 



LOoEp 

nW 



V 



/So 
-w -v/T 



d^ 



1 



w 



VI - {^/wy {Ek + Wo - 0^ 
d^ 1 



ko VI - i^/wr i^k - wo - 0' 

For the incoherent part of the spectral function we find 



A''=(fc,w < -Wo) = 



and 



1 wogp [W^ - (w + wo + fif] ^ J% d^(5(w + wo + M - 
TT [W^ -(aj + uJo + Ai)2] {oj + n-^k-'R.e S(fc, uj)f + {uJoEp)- 



(12) 



(13) 



1 woi5;p[PF2_(^_^^ + ^,)2]. 

A'^(fc,w>wo) = 5 . 14 

7^[W^-(cu-u;o + ^ir]{LU + |2-^k-Rcl:{k,u;)f + {u;oEp)^ 

Finally, the chemical potential iJ for a given electron density n is determined by 



1 /""^ 

AtX] / dw^(A;,w)nF(w) = n. 



(15) 



2.2.2. Strong coupling Hamiltonian (1) with the coefficients (3) represents the 
Hamiltonian of small polarons, which are the correct quasiparticles in the SC limit. 
Using the procedure outlined above, we obtain the polaron self-energy as 



W 



w 



+ 



l-nF(^ + 7?) 
w - swo - i^ + v) ' w -I- swo - i^ + v) 



(16) 
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ReT,{k,u;) 
and 

ImE(A:, w) 



--f{k,^,s) 



+ 



<-v) 



"9" 2^ 



u) - su>o - i^ + v) ^ + stoo - i^ + v) 



(17) 



S>1 



+9{-uj - suJo)d[uj + su>o - i^ + v)]} ■ 
Here ^(a;) is the Heaviside step function and we have used the definition 



f{k, ^, s) {e{uj - sujo)S[oj - s(Jo - K + v)] 

(18) 



2(X 
W 



1 



2\s 



S\ S\ 

The coherent part of the spectrum, non-zero for < wq, is given by 

A%k,u;) = ZkS[u; - {Ek + v)] 
with the renormaHzed band energy Ek being the solution of 

Ek^Ck + Rc^{k,Ek+v), 
where = —W cosk, and the spectral weight takes the form 



cos(2fc) . (19) 



(20) 
(21) 



Zu = 



_{Ek-sujo-0^ ' {Ek + suo-0^\ ■ ^^^^ 
The imaginary part of the self-energy, determining the incoherent excitations, is 
non-zero only for |a;| > wq. Consequently, we get 



Iml](fc,w ^ T'^o) = 
with 



W 



^9{t<^ - sujq) 



s>l 



f^{k,^,s) 



w 



w 



d^6{uj ± swo - r? - > (23) 



1 



to ± SLOq — 77 

w 



(24) 



If ojo > 2W, for a given w, only one term of the sum in equation (23) contributes. The 
corresponding index a is determined by the conditions 



or 



_ d^5{L0 — (TUjQ — 77 — ^) = 1 
-w 

6{—u) — auo) = 1 and / _ d^S{LO + auo — rj — ^) = 1. 

J-w 



(25) 



(26) 
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2lT 



(27) 



Here we have defined 

f'^{k,w,s) = f{k,^ = uiTswo-v,s). (28) 
The polaron spectral function determines the equilibrium properties of the 
spinless HM in the SC regime. The photoemission spectra, however, are 
determined by the electron spectral function which is related to the retarded Green 
function containing electronic operators. According to the canonical Lang-Firsov 
transformation [20], which defines small polaron states, the relation between the 
polaronic operators Ci entering the Hamiltonian (1) with the coefBcients (3), and 
the transformed electron operators Ci, reads 

= cxp[,g(6] - b^)]c^ , = exp[-9(6t - 6.)]c] (29) 
or, in the Bloch representation. 



(30) 



J 3 

We start our derivation of the relation between the polaronic and electronic 
spectra from the time-ordered Green function [26] for the (transformed) electron 
operators 

(31) 



and factorizc the statistical averages with respect to polaron and phonon variables 
G^{kMM) = ^e'^'^GT(t,d)(Tte^[''^(')-''.('51e-'^l'''+^(°)-''.+^(°)]) , (32) 

d 

where t = ti — t2 and G'^{d, t) = —i{Tic- (t)ct^j^(O)) represents the time-ordered Green 
function of polaron operators fulfilling 

GT(d,t) =iV-i^e-"='''GT(fc',t), G^{k',t) = -i{Ttc^,{t)cl,{Q)). (33) 

k' 

The averages over the phonon variables in equation (32) will be (;valuated using 
mutually independent local Einstein oscillators having the time-dependence hj{t) = 
Q-iujotjj. 'Working in the low-temperature approximation we obtain 



GT(fc,w) =e-» G'^{k,uj)+e-3 



S>1 



s! N 



(34) 



{u}-\-8U}Q—ie)t 



/ dtG>{k',t)e'^'^-"^°+'''>^ + / dtG<{k',t)e' 

f., Uo J-oo 

with G-^{k,t) = — i(cj.(t)c^(0)) and G'^(k,t) = i{cl.{0)Cf.{t)) , and the convergence 
factor exp(— e|t|), e — > 0+. Introducing the generalized function ({iv) = [ui + ie]~^, 
e 0+, we get 



GT(fc, w) = e-B"G^{k, Lo) + e-s' ^ 



E 



s\ N 
>1 h 



du'G>{k',u;'^ 



/OO 
-OO 



(35) 
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The Green functions G^{k,uj) are related to the polaron spectral function Ap{k,uj) 
through 

G<{k,Lu) =inF{Lo)Ap{k,Lo), G>{k,io) = - nF{uj)]Ap{k,u) . (36) 
At T = 0, we have 



G'^{k,w)=e-<''G'^{k,Lj)+e-o'Y^ 



E 



s! N 

>1 k 



dijo'Ap{k',ui')({co — suio — w') 



+ / duj'Ap{k',J)C{oj + soJo-uj') 

J —OO 



(37) 



The relation between the time-ordered Green function G'^{k,co) and the associated 
retarded Green function G^(fc,w) in the low-temperature approximation reads 

M 

and hence we obtain 



ImG^(fc,w) = -^ImG'^{k,Lo) 



(38) 



Ap{k,u>) = --■^ImG'^(fc,a;) . 

TT \UJ\ 



(39) 



Of course, equations (38) and (39) hold for the electron Green function as well. 
Consequently, the electron spectral function Ae{k,co) is expressed in terms of the 
polaron spectral function Ap{k,Lj) as 

Ae(fc, w) = e-^'Ap{k, w) + e-«' ^ Yl 



si 



Y^[Ap{k' ,uj — sujo)(^{<^ — sujo) + Ap{k' ,00 + sLUo)9{—oj — sluq)] . (40) 



Similar results have been derived before in [27, 28]. 



2.2.3. Intermediate coupling As we shall see in section 3, the results of the WC (SC) 
approximation are in good agreement with numerical calculations [18, 19] if A, g <C 1 
{\,g ^ 1). However, the cross-over between these limiting cases, revealed by the 
numerical calculations, appears to be out of reach for the analytical formulae hitherto 
deduced. 

To interpolate between WC and SC, we shall modify the method of canonical 
transformation by Lang and Firsov [20]. In the latter, the term of the HM (1) linear 
in the local oscillator coordinate Xi = xo{bl + b^) is completely eliminated by the 
translational transformation U = exp[^^ gicjc^ (6| — b^)]. As a result, the local lattice 
oscillator at the site i is shifted by Aa;^ = 2a;o5cJq, if occupied by a charge carrier, 
whereas there is no such deformation at unoccupied sites. To generalize this picture, 
we abandon the site localization of both the charge c;arrier and the lattice deformation 
in the transformation. Physically, these localizations will be destroyed with increasing 
hopping rate and charge-carrier concentration in the IC regime. Different canonical 
transformations, taking into account charge-density-wavc order at n = 0.5, have been 
proposed, e.g., by Zheng et al [29]. However, in this approach, there is no filling 
dependence of their variational parameter and of the mean lattice deformation, which 
is crucial for a correct description of the adiabatic IC regime. 

The probability for the charge carrier to be found at a distance \rj — rj| from 
the center r, of the polaron will be assumed to be proportional to exp(— |rj — r^l/i?). 
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In one dimension, and setting the lattice constant to unity, we use the normahzed 
distribution 

=e-^tanh^. (41) 

Accordingly, the shift of the local oscillator with coordinate Xi is assumed to be 

Axi = 2xQgKclci + 2xQgj{l - c\ci ) (42) 

with 

ii: = tanh;^, ^ = 2Ke~Tin. (43) 
2R 

The last term in equation (42), characterizing the mean lattice deformation 
background, takes into account the influence of nearest-neighbor sites only. 
The canonical transformation leading to the oscillator shift (42) reads 



U = exp 



7 = - 7 . (44) 



Carrying out the transformation U^HU = H for the Hamiltonian (1) with (2), the 
terms of H containing polaron operators are modified with the following coefficients 
[cf. equations (2) and (3)] 

Cu = gujoi^ - 7)(&I + q,,) = te-^^^''-'^ -'1+'^^ . (45) 

Moreover, we have r] = —jj, — Ep[y(2 — 7) + 27(1 — 7)]. 

Owing to numerical problems which occur in certain parameter regimes when 
using the full Hamiltonian with the coefficients (45), the variational parameter R of 
the transformation will be determined in the first approximation, which is analogous 
to the Hartree approximation. The corresponding polaron ^ectral function is given 
as Ap{k,Lj) = 6[uj — (^fe + 77)], where = —Wcosk with W = W exp[— (751)^], and 
T] as defined above. R is then defined by the position of the minimum of the total 
energy per site E/N in the first approximation, i.e.. 



= I 0{-^ - r?) + /xn + Epj^ , (46) 

^ ttW J-w _ (^^/wy 



with the condition for /j. 



' ^"de4i=i=i-. (47) 



ttW J-w ^i_(^/H^)2 



The last term in equation (46) arises in H from the lattice deformation background 
at finite concentration n. 

R determined in this way for each set of model parameters will be used to calculate 
both the electron and polaron spectral function, taking into account the multi-phonon 
processes included in H [cf. equation (45)]. 
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Using the same procedure as in the SC case, we calculate the self-energy and 
spectral function at T = 0, finding 



ReS(A;,w) 



W 
2^ 



S>1 



0{^ + v) 



f{k,^,s) 



LJ - suo - i^ + ri) + '^^0 - + ??) 



-27(1-7)^7' 

TT 



W 



dC 



0{i + ri) 



cosk 



uj - uja ~ {S, + ri) Lu + luq - + Tj) 



-\2 EpUJo , 



ttW J-w 



/w 
-w 



dC 



1 - {i/W? 



+ 



Wo - + w + Wo - + ^) 



(48) 



and 



W 



ImE(fc,a; ^ T^o) = -^^0{T(^ - sujq) '' /JlJ / „ d^S{uj ± s(Jo - - 



W 



+ cosfc 



T27(l-7)i^p(X±)- 

/w 
_ d^(5(a; ± wo - - ^) 



J-w 



(49) 



Here W = Vl^e"^'''^) and /, / are defined as in equations (19) and (28) but with g 
replaced by 57. Note that the above equations reproduce the corresponding SC results 
in the limit = (7 = 1), and also the WC ones in the limit R = 00 (7 = 0). 

The coherent spectrum is again given by equations (20) and (21), with the spectral 
weight 



z^' = |1 - [dReJ:{k,uj)/duj]^=E,+^\ , 
whereas the incoherent part takes the familiar form 

1 ImE(A:,a;) 



A-(fc,c.) 



(50) 
(51) 



7t[lo- (a RcE(fc, c^)]2 + [ImE(fc, lu)]^ ' 

Finally, equation (40), which determines the relation between the electronic and 
polaronic spectral functions, also applies to the present case if g is replaced by 57 
throughout. 



3. Numerical results 



As in section 2.2, we first discuss the WC and SC limits, before turning to the 
important IC regime. Since we use a finite number of momenta k, it is not possible 
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Figure 1. (colour online) Coherent {A'^, ) and incoherent {A"', ) parts 

of the spectral function in the weak-coupling approximation for different band 
fillings 71. Here ujo/t = 0.4 and Ep/t = 0.1. 



to tune the band filling n (via the chemical potential /i) to a specific, desired value 
with arbitrary accuracy. In order to simplify the discussion of the different density 
regimes, we therefore report rounded values of n in the figures and the text. The 
largest deviations of the actual n from the value reported occur in the SC case for 
which, however, the density dependence is very weak (see section 3.2). 

3.1. Weak coupling 

Figures 1-3 show the electronic spectral function A{k,uj) obtained from the WC 
approximation. The coherent spectrum (ImS = 0) is given as the solution of 
equations (10)-(12) with energies \Ek — /i| < ujq. For \uj\ > ujq, A^'^, calculated 
according to equations (13) and (14), consists of peaks having widths proportional 
to Ep. A comparison of figures 1 (a) and 2(a) shows the spreading of the coherent 
spectrum with increasing luq. Finally, comparing figure 1(d) [2(b)] with figure 3(a) 
[3(b)], we observe a broadening of the peaks in with increasing Ep. 

The evolution of the spectral function with increasing carrier density n is 
illustrated in figure l(a)-(d). The coherent part is shifted inside the spectrum as 
a function of n (i.e., with n). Additionally, the shape of A"^ is also affected by n, due 
to the dependence of equations (13) and (14) on the chemical potential fi. 

In figure 4, we plot the total spectral weight J dk J dujA{k,uj), as well as the 
coherent weight J dk J dujA'^{k,uj). Note that for any Ep > 0, the coherent band Ek 
is restricted to the interval [— a;o,^^o]i so that the corresponding coherent weight is 
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Figure 3. (colour online) As in figure 1, but for Ep/t = 0.5 and fixed band filling 
n = 0.4. Here (a) too/t = 0.4 and (b) uio/t = 2. 



significantly smaller than the Ep = value of unity. This reduction is less pronounced 
for larger phonon frequency ujo/t — '2 ~ W/t (figure 4). Furthermore, we see that with 
increasing Ep, the sum rule for A(k,u!) becomes more and more violated, as expected 
for a WC approximation (for a more detailed discussion of sum rules see section 3.2). 

Finally, figure 5 displays the coherent band dispersion Ek at small carrier density 
n — 0.1. As known from the single-polaron problem, the coherent weight Zk drops to 
zero as the bare phonon dispersion intersects with the renormalized band. This gives 
rise to a flattening of the coherent band at large k [30], which is well reproduced by 
the simple WC approximation. 

3.2. Strong coupling 

We now turn to the opposite, SC limit. The theory presented in section 2.2.2 directly 
yields the polaronic spectrum Ap(k,uj), results for which are shown in figure 6(a) for 
ujo/t = 0.4 and Ep/t = 4. The spectrum is dominated by a coherent polaronic band 
with negligible width (for the dependence of the bandwidth on Ep see figure 7) having 
a spectral weight Zk close to unity (cf. figure 8). This suggests that small polarons 
are the correct quasiparticles in the SC regime. Note that opposite to the WC case, 
where the sum rule for the spectral function becomes more and more violated with 
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Figure 4. Coherent ( ) and total { ) spectral weight from A{k,u>) (see 

text) for u>o/t = 0.4 (□, I) and wo/t = 2 (O , • ) as a function of the polaron 
binding energy Ep. Here ^ = 0. 








Figure 5. Renormalized band for uJo/t = 0.4, Ep/t = 0.1 and n = 0.1. 

Symbols indicate non-zero coherent spectral weight Zk ■ The horizontal line ( ) 

corresponds to {Eq +uJo)/t. 




Figure 6. (colour online) Coherent {Ap, ) and incoherent {Ap, ) 

parts of the polaron spectral function in the strong-coupling approximation. Here 
n = 0.4, Ep/t = 4 and (a) uJo/t = 0.4, (b) uJo/t = 2. 
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-□0000000^ 
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Figure 7. Renormalized band for n = 0.4, u>o/t = 0.4 and different values 
Ep/t = 2 (O ) and Ep/t = 4 (□). 




Figure 8. Coherent (polaronic, ) and total ( ) spectral weight from 

Ap{k,u>) for LUo/t = 0.4 {□, I) and uq/* = 2 (O , • ), as a function of the 
polaron binding energy Ep. Here /i = 0. 




Figure 9. (colour online) As in figure 6, but showing the electronic spectral 
functions ( ) and A]^ ( ). 
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Figure 11. (colour online) As in figure 10, but for ujo/t = 2. 



increasing coupling (figure 4) , here the SC approximation becomes increasingly better 
with increasing Ep (figure 8). 

We would like to point out that such changes in the total spectral weight are 
absent in the work of Alexandrov and Ranninger [27], since the latter was restricted 
to the lowest (first) order of the self-energy, similar to the Hartree approximation 
discussed in section 2.2.3. In general, the total spectral weight contained in the 
(electronic or polaronic) spectral function for given parameters depends on the 
approximations made. As illustrated by figures 4 and 8, the total spectral weight 
approaches the exact value of unity in the WC and SC regimes, respectively, so 
that the normalization of the spectrum serves as a measure of the validity of the 
underlying approximations. Since in the present case even the first moment (i.e., the 
normalization) shows deviations from exact results, we have refrained from checking 
the more complicated sum rules derived in [31]. This is also true of the IC case 
discussed below. 

The effect of increasing the phonon frequency loq can be seen by comparing 
figures 6(a) and (b). Most noticeably, for larger ojq, the width of the coherent polaron 
band — roughly scaling proportional to — is larger. 

The polaronic spectrum is related to the electronic spectrum Ac by equation (40), 
and typical results in the adiabatic and non-adiabatic regimes are shown in 
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Figure 12. Total energy per site E/N as a function of the variational parameter 

R for different values of Ep from the Hartree approximation. Here n = 0.4 and 
(a) ujo/t = 0.4, (b) ojo/t = 2. The values of Ep/t are (from top) (a) 0.1, 1, 2, 2.5, 
3, 4, (b) 0.1, 1, 2, 4, 6, 8, and 12. 



figure 9. Although strictly speaking a distinction between coherent and incoherent 
contributions cannot be made in the case of [cf. equation (40)], it is useful to 
separate the two terms on the r. h. s. of equation (40), and to identify the first as the 
contribution of the coherent polaron band. 

Carrying out the transformation from Ap to Ae according to equation (40), the 
weight of the coherent polaron band visible in figure 6 approximately acquires a 
prefactor e"^ . The remaining contributions to the electronic spectrum correspond to 
phonon-assisted photoemission processes. 

In the case of A^, the main difFercnce between the adiabatic [figure 9(a)] and the 
non-adiabatic regime [figure 9(b)] is the significantly larger weight of the coherent band 
for large lvq since g"^ = Ep/uJo. Consequently, the weight contained in the incoherent 
excitations is noticeably reduced. 

The results in figures 10 and 11 show a certain dependence on the band filling n, 
but there occur no qualitative changes even at IC Ep/t = 2. This is in contrast 
to recent numerical work [18,19]. In particular, the spectrum in figure 10(b) 
is substantially different from figures 14(c) and (d), which are all for the same 
parameters. We shall see below that a more satisfactory description of the real physics 
can be obtained using the variational approach discussed in section 2.2.3. 

3.3. Cross- over from weak to strong coupling 

As seen in the preceding sections, both the WC and SC approximations are not 
capable of accounting for the recently discussed carrier density-driven cross-over from 
a polaronic system to a metallic system with phonon-dressed electrons [18, 19]. In fact, 
the electronic spectrum always remains WC/SC-like in character. In order to obtain 
a reasonable, analytical description of the IC regime, we therefore use the variational 
approach proposed in section 2.2.3. 

Figure 12 shows results for the total energy per site as a function of the variational 
parameter R for different values of Ep. The evolution is very similar to the large-to- 
small polaron cross-over in the one-electron case [25]. Note that the present results 
have been obtained using the Hartree approximation, as discussed in section 2.2.3. 

For weak coupling Ep, in the adiabatic case [figure 12(a)], we find a minimum in 
the total energy at a finite value of R. Upon increasing Ep, a second local minimum 



Spectral functions of the spinless Holstein model 



16 




Figure 13. (a) Total energy per site E/N as a function of Ep obtained from 
A{k,u!) (see text), where n = 0.4 and u>o/t = 0.4. Also shown are quantum 
Monte Carlo results at inverse temperature /3t = 8 obtained with the method of 
[18]. (b) Exact diagonalization (ED) and density matrix renormalization group 
(DMRG) results for the renormalized kinetic energy as a function of carrier density 
n for Ep/t = 2, ujo/t = 0.4 and different cluster sizes N. 




Figure 14. (colour online) Spectral functions from the variational approach 
for different values of Ep, n = 0.4 and ujo/t = 0.4. The parameters are (a) 
Ep/t = 0.1, R = 1.5, (b) Ep/t = A, R = 0.1 and (c) Ep/t = 2, R = 1.3. The 
values of R have been chosen according to figure 12(a). Panel (d) shows cluster 
perturbation theory results for photoemission [A~ (k,u!), dashed lines] and inverse 
photoemission [A~^ {k,ijj)] for n = 0.4 and -Ep = 2 (taken from [19]). 




Figure 15. (colour online) Spectral functions from the variational approach for 
different values of Ep, n = 0.4 and ujo/t = 2. The parameters are (a) Ep/t = 0.1, 
R = 0.7, (b) Ep/t = 8, R = 0.2, (c) and (d) Ep/t = 4, R = 0.5. The values of R 
have been chosen according to figure 12(b). In the white regions in (d) < 
(see text). 



starts to develop near R — 0, associated with the small-polaron state which becomes 
the ground state in the SC Umit. The jump of the optimal value of R from a large to 
a small value at a critical Ep suggests a first-order transition from an extended to a 
small-polaron state. However, such a sharp transition is absent in the single-polaron 
case, and not expected for n > either. The discontinuous cross-over appearing in 
our results for wq/^ ^ 1 is a consequence of the approximation used. 

In contrast, in the non-adiabatic regime [figure 12(b)], there exists only a single 
minimum, which shifts to smaller R with increasing coupling Ep. Moreover, compared 
to figure 12(a), the dependence of the total energy on R is much weaker. 

Apart from the comparison of the spectral function with other methods presented 
below, variational approaches — often not capable of yielding dynamic properties — are 
usually judged by the total energy as opposed to exact data. Figure 13(a) shows 
the second-order results for the total energy as a function of Ep, using the optimal 
values of the parameter R as determined from the Hartree approximation. Clearly, 
the agreement between our IC approach and results from quantum Monte Carlo [18] 
is very good at weak and strong coupling (the variational approach reproduces the 
WC and SC limits of sections 2.2.1 and 2.2.2), whereas there are notable deviations 
at IC. The missing higher-order corrections — causing the violation of the sum rule 
discussed earlier — lead to generally overestimated values of the energy. Note that the 
agreement of the energy with exact results is even better in the non-adiabatic regime 
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(not shown). 

To illustrate the density-driven cross-over from (large) polarons to slightly dressed 
electrons (scattered by diffusive phonons), figure 13(b) reports exact numerical results 
for the renormalized kinetic energy. The latter may serve as a measure for the 
carrier mobility. The increase as a function of n may be interpreted as originating 
from the overlap of the displacement clouds surrounding the carriers and, finally, the 
dissociation of the polaronic quasiparticles at large n. 

We again begin the discussion of the spectral functions with the adiabatic case 
LOo/t = 0.4. Following the discussion of section 2.2.3, we determine the optimal R 
from the position of the minimum of the total energy for given Ep and n (figure 12). 

Figures 14(a) and (b) resemble closely to the WC and SC results of figures 1(d) 
and 10(b), respectively. However, for IC Ep/t = 2 [figure 14(c)], we find a rather 
metallic spectrum with a broad main band crossing the Fermi level, and with low- 
energy excitations available. The corresponding polaronic spectrum (not shown) 
reveals that the variational approach correctly predicts the absence of well-defined 
polaronic quasiparticles, as suggested by the non-negligible incoherent contributions 
in Ap{k,Lu), lying close in energy to the coherent band. This is in contrast to the 
SC case, where polaronic quasiparticles dominate, and in which the coherent band is 
well separated from the incoherent excitations. Note that in figure 14(c) [and also 
in figures 15(c) and (d)] the incoherent part becomes slightly negative for uj > —W 
{lu < W) at large (small) k, which is an artifact of our approximation. Remarkably, the 
overall features of the spectrum in figure 14(c) are very similar to the corresponding 
numerical results in [19], reproduced in figure 14(d). 

Finally, in the non-adiabatic regime, it has been found numerically that small 
polarons remain the correct quasiparticles even at large band fillings n [17, 18]. Again, 
the variational approach is able to describe the physics correctly. In particular, 
figures 15(c) and (d) reveal that a dominant coherent band, well separated from 
incoherent excitations and having a relatively small bandwidth, exists even at IC. 

As for the WC and SC cases discussed above, we have also checked the total 
spectral weight in the present case. We find that the deviations from the exact value of 
unity are largest for WC, whereas the sum rule is fulfilled to within 1-10% (depending 
on ojQ and Ep) at IC and SC. 

4. Conclusion 

We have presented an analytical treatment of the one-dimensional spinless Holstein 
model based on calculations of the self-energy in the framework of the generalized 
Matsubara functions. To connect the analytical results to previous numerical ones, 
the electronic spectral function determining the photoemission spectrum has been 
computed for finite carrier concentrations in dependence on the electron-phonon 
coupling strength and the phonon frequency. 

In the strong-coupling limit, the electron spectral function has been deduced 
from the spectral function of small polarons. However, it was shown that the electron 
picture and the small-polaron picture both cease to be correct if we approach the 
intermediate-coupling regime from the weak- and strong-coupling side, respectively. 

To describe the cross-over from the strong- to the weak-coupling limit, 
an interpolation scheme based on a generalization of the Lang-Firsov canonical 
transformation has been proposed. The latter was defined for each set of model 
parameters by a distance R which characterizes the charge distribution across the 
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polaron volume. The parameter R, as deduced from the minimum of the total energy 
in the first approximation, was shown to increase with decreasing Ep corresponding 
to a cross-over from large to small polarons — and, at the same time, the coupling 
dependence of R was found to be stronger in the adiabatic than in the non-adiabatic 
case. 

The spectral functions calculated at weak, strong and intermediate coupling are 
in a good agreement with recent numerical calculations. Moreover, analytical results 
deduced by means of the self-energy calculations enabled us to distinguish between the 
coherent and incoherent parts of the spectrum. Most importantly, starting from the 
strong-coupling limit, it was shown that the spectral weight of the incoherent polaron 
spectrum increases with decreasing coupling Ep, and that the energy separation of the 
incoherent peaks from the coherent spectrum is continuously reduced in the adiabatic 
regime. On the contrary, the coherent part of the electronic (photoemission) spectra 
is reduced by a factor e~* , and the incoherent part representing the phonon-assisted 
photoemission processes becomes increasingly dominant with increasing coupling. 

Finally, at intermediate coupling and finite carrier densities, our results support 
recent numerical findings which suggest that the system can no longer be described 
in terms of (small) polaronic quasiparticles [18, 19]. 
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